Extending the range of the inductionless magnetorotational instabihty 
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The magnetorotational instability (MRI) can destabilize hydrodynamically stable rotational flows, 
thereby allowing angular momentum transport in accretion disks. A notorious problem for MRI is 
its questionable applicability in regions with low magnetic Prandtl number, as they are typical for 
protoplanetary disks and the outer parts of accretion disks around black holes. Using the WKB 
method, we extend the range of applicability of MRI by showing that the inductionless versions 
of MRI, such as the helical MRI and the azimuthal MRI, can easily destabilize Keplerian profiles 
oc r~^^^ if the radial profile of the azimuthal magnetic field is only slightly modified from the current- 
free profile (X . This way we further show how the formerly known lower Liu limit of the critical 
Rossby number, Ro ~ —0.828, connects naturally with the upper Liu limit, Ro ~ +4.828. 
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Initiated by the seminal work of Balbus and Haw- 
ley [l[ , the magnetorotational instability has become the 
standard explanation for turbulence and enhanced angu- 
lar momentum transport in accretion disks around black 
holes and proto-stars. While MRI is thought to be a ro- 
bust phenomenon in the hot parts of accretion disks, a 
notorious problem concerns the viability of MRI in other 
regions, such as the outer parts of black hole accretion 
disks y and the "dead zones" of protoplanetary disks 
Q. This has to do with the fact that the onset of MRI 
demands that both the rotation period and the Alfven 
crossing time in vertical direction are shorter than the 
timescale for magnetic diffusion For the case of a 
vertical magnetic field applied to a disk of height 
H this means that both the magnetic Reynolds number 
Rm = ^QaH^n and the Lundquist number S = ^quHva 
must be larger than one, and that S < Rm (fi is the 
angular velocity, /iq is the magnetic permeability, a the 
conductivity, va '■— Bz/ s/Jm)P is the Alfven velocity, 
with p denoting the density). In a disk with given 
size, angular velocity, and magnetic field strength it is 
then often the spatially varying magnetic Prandtl num- 
ber Pm = v/r], i.e. the ratio of viscosity v to magnetic 
diffusivity 77 := (/xqct)"^, that determines the values of 
Rm and S, and hence the fate of MRI. 

For the case without external Bz things are even more 
complicated since the MRI-triggering magnetic field, in 
this case dominated by the azimuthal component B^, 
must be produced in the disk itself, very likely by some 
sort of an a—^l dynamo process [5(]. This combined, loop- 
like action of MRI and self-excitation has attracted much 
attention in the past, with many open questions concern- 
ing issues of numerical convergence [6] , as well as the role 
of disk stratification and vertical boundary conditions 
. Again, the most interesting case appears in the limit 
of low Pm. While Lesur and Longaretti ^ have argued 
for a power-law decline of the turbulent transport with 
decreasing Pm, there are also indications for the exis- 



tence of some critical Rm in the order of 10'^...10'* for the 
MRI-dynamo loop to work [lo| . 

Exactly this situation, characterized by low Pm and 
a significant or even dominant B^, is the subject of in- 
tense theoretical and experimental research initiated by 
Hollerbach and Riidiger For the ratio of to 

Bz being on the order of 1 and B^{r) cx l/r, helical 
MRI (HMRI) was shown to work also in the induction- 
less limit [12], Pm = 0, and to be governed by the 
Reynolds number Re = RmPm^^ and the Hartmann 



number Ha — SVi 
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quite in contrast to standard 



MRI (SMRI) that is governed by Rm and S. 

Somewhat disappointingly, a crucial limitation of this 
surprising kind of MRI was identified by Liu et al. [13] 
who used a WKB approach to find a minimum steepness 
of the rotation profile, expressed by the Rossby number 
Ro := r{2n)~^dn/dr < Rolll = 2(l-\/2) « -0.828. 
This limit, which we call lower Liu limit (LLL) in the 
following, implies that the inductionless HMRI in the 
case when B^(r) cx l/r does not extend to the most 
relevant Keplerian case, characterized by Roxep = —3/4. 
In addition to the LLL, the authors found also a second 
threshold of the Rossby number, which we call the upper 
Liu limit (ULL), at Roull = 2(1+^2) « +4.828. This 
second limit, which implies a magnetic destabilization of 
extremely stable flows with strongly increasing angular 
frequency, has attained nearly no attention up to present, 
but will play an important role below. 

The existence of the LLL, together with a variety of 
further predicted parameter dependencies, was confirmed 
in the PROMISE experiment working with a low Pm liq- 
uid metal [l^ . Present experimental work at the same 
device aims at the characterization of the azimuthal MRI 
(AMRI), a non-axisymmetric "relative" of the axisym- 
metric HMRI, which is expected to dominate at large 
ratios of Ba to Bz la. However, AMRI as well as in- 



ductionless MRI modes with any azimuthal wavenumber 
(which may be relevant at small values of Bff,/Bz), seem 
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also to be constrained by the LLL as recently shown in 
a unified WKB treatment of all inductionless versions of 
MRI [ill- Actually, it is the apparent failure of HMRI, 
and AMRI, to apply to Keplerian profiles that has pre- 
vented a wider acceptance of those inductionless forms of 
MRI in the astrophysical community. Only recently, the 
intricate, though continuous, transition between SMRI 
and HMRI was explained in some detail by showing that 
it involves a spectral exceptional point at which the in- 
ertial wave branch coalesces with the branch of the slow 
magneto-Coriolis wave • 

Given the fundamental importance of whether any sort 
of inductionless MRI could possibly work in the low Pm 
regions of accretion disks, it is quite natural to ask for 
how to extend the range of its applicability beyond the 
LLL. In a first attempt, the stringency of the LLL for 
BAr) (X 1/r was questioned by Riidiger and Hollerbach 
[18| who had found an extension of the LLL to Keple- 
rian values in global simulations when at least one of the 
radial boundary conditions was assumed electrically con- 
ducting. Later, though, by distinguishing between con- 
vective and absolute instabilities for the travelling waves 
such as HMRI, the LLL was vindicated even for such 
modified electrical boundary conditions [T^]. A second 
attempt was made in 19] treating HMRI for non-zero, 
but low S. It was found that for Btjj{r) cx l/r, the es- 
sential HMRI mode extends from S* = only to a value 
S « 0.618, and allows for a maximum Rossby number of 
Ro ~ —0.802 which is indeed slightly above the LLL, yet 
below the Keplerian value. Close to this critical point, 
the essential HMRI is then replaced by a helically mod- 
ified SMRI. A third possibility arises by noting that the 
saturation of MRI could lead to modified flow structures 
with parts of steeper shear, sandwiched with parts of 
shallower shear [20| . 

In this Letter, we discuss another promising way of ex- 
tending the range of applicability of the inductionless ver- 
sions of MRI to Keplerian profiles, and beyond. Rather 
than relying on modified electrical boundary conditions, 
or on locally steepened ri(r) profiles, we will evaluate 
B^{r) profiles that are shallower than l/r. The main 
idea behind that is the following: Assume that in a low- 
Pm region, characterized by S* << 1 so that standard 
MRI is reliably suppressed, Rm may still be sufficiently 
large for inducing azimuthal magnetic fields, either from 
a prevalent axial field Bz or by means of a dynamo pro- 
cess without any pre-given Bz- If i?0 is produced exclu- 
sively by an isolated axial current, we get oc l/r. The 
other extreme case, B^ cx r, corresponds to the case of 
a homogeneous axial current density in the fiuid which 
is already prone to the kink- type Tayler instability plj . 
even at Re = 0. For real accretion disks with compli- 
cated conductivity distributions in radial and axial di- 
rection, quite a variety of intermediate B^{r) dependen- 
cies between oc l/r and oc r profiles is well conceivable. 
Leaving those details aside, here we focus on the generic 



question which deviations of the B(j,{r) profile from l/r 
could make HMRI (or AMRI) a viable mechanism for 
destabilizing Keplerian rotation profiles. By defining an 
appropriate magnetic Rossby number Rb we will show 
that the instability extends well beyond the LLL, even 
reaching Ro = when going to Rb — —0.5. Evidently, 
in this extreme case of uniform rotation the only avail- 
able energy source of the instability is the magnetic field. 
Most interestingly, by tracing the instability threshold 
further into the region of positive Ro in the Ro — Rb 
plane, we find a natural connection with the ULL whose 
meaning was a somewhat mysterious conundrum up to 
present. 

We set out from the equations of incompressible, 
viscous and resistive magnetohydrodynamics, i.e. the 
Navier-Stokes equation for the velocity field u and the in- 
duction equation for the magnetic field B, together with 
the continuity equation for incompressible flows and the 
divergence-free condition for the magnetic field: 
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B- Vu-u - VB + T^V^B, (2) 
0, V • B 0. (3) 



We consider a purely rotational flow exposed to a mag- 
netic field comprising a constant axial component and an 
azimuthal one with arbitrary radial dependence: 

Uo(r)=rr!(r)e^, Bo(r) = i?0(r)e^ + S°e,. (4) 

To study flow and magnetic field perturbations on this 
background we linearize the equations in the vicinity of 
the stationary solution by assuming u = Uq + u' , p = 
Pq + p', and B = Bq -I- B' and leaving only terms of first 
order with respect to the primed quantities. Introducing 
the total wavenumber |kp — k"^ + k^, and a — /cz/lkj, 
where fc^ and kz are the radial and axial wavenumbers 
of the perturbation, we define the viscous, resistive, and 
two Alfven frequencies corresponding to Bz and B^: 



j^=77|k|- 
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(5) 



Then, we define the ratio (3 of the two field components, a 
re-scaled azimuthal wavenumber n, the Reynolds number 
Re, and the Hartmann number Ha as follows: 

p — a n = — , Re — a — , Ha = ——^ . [b) 

The steepness of 5^(r) will be measured by the hydrody- 
namic Rossby number, and the steepness of B^(r) by the 
corresponding magnetic Rossby number: 



r do. r duJA^ 

^""-m^' ^^-2u^A. dr 



(7) 
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By employing the same short-wavelength (WKB) ap- 
2^ . but now including Rb, we end 
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proximation as m 
up with a system of 4 coupled equations for the pertur- 
bations of arbitrary azimuthal wavenumber, yielding the 
ultimate dispersion relation det (M — A/) — 0, with A 
denoting the (complex) growth rate in units of aD, and 
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where Rm = RcPm is the magnetic Reynolds number. 
As a first test case, this relation can be applied to the 
kink-type Tayler instability that has recently been ob- 
served in a liquid metal experiment j2]J . In the relevant 
limit with Pm = and Re = we deduce from the Bil- 
harz criterion [23i] the following condition for marginal 
stability: 

^ (l + Ha^M + l)^)^-4HaVM + l)2 

4Ha2/?2(l +Ha2(n/3 + 1)2) ' ^ ' 

For Rb = 0, which corresponds to cx r, and taking 
the limit /3 — )• oo, we obtain /3Ha = (1 - (1 ± n)2)~i/2 
which would become equal to 1 for n ~ =pl. Trans- 
lated to the real experiment with « 2.4/r and a very 
rough estimate kr ~ t/t, we find a value of Haoxp ■— 
B(f,(r)ryJ <y I pv « 34 which is not too far from the exper- 
imentally observed value of 22 [2l[ . 

Our main focus here is, however, on the limit Re — >■ oo 
and Ha ^ oo that is relevant for MRI. Assuming for the 
moment Pm = (which will be slightly relaxed later), 
and inserting the optimal relation between Re and Ha, 



Re = 2RbV3Rb + 2(\/l + 2Rb + (10) 

(obtained in the manner described in [l6[), we find from 
the Bilharz criterion ^231 the dependence of the critical 
Rossby number on Rb, n, and /3: 



Ro± = -2 



2/32(n/3-M)2 



(11) 



where F = (n^ -I- 1)^ - 2/32Rb. Note that under the 
assumption Pm — the dispersion relation possesses 
an exact solution, which after being expanded into the 
Taylor series with respect to the interaction parameter 
N = Ha^Re~^ in the vicinity of N = 0, is 



A 



(n ± 2Vl + Ro) - Re"^ 
^(Ro-F2)(n;3-f 1) 



N F± 



\/l + Ro 



+ 0(N2). (12) 



At n = 0, Rb = — 1 and Re — >■ c» the growth rates 
reduce to those derived in [l^. In the limit N — > and 



Re — > oo the stability boundary is obtained when the 
real part of the term linear in N vanishes. This condition 
leads exactly to equation (fTTj) . which also confirms the 
correct application of the Bilharz criterion. 

With the goal to find extremal values of Ro that are 
compatible with marginal stability, we can further opti- 
mize /3 and n (or ot) according to 



/3opt — 



-1 



n±y^2Rb' 



^opt 



) ^ (13) 
) y=2Rb 



to obtain Ro^pt(Rb)= - 2-4Rb±2(2Rb(2Rb+l))i/2^ or 
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1 (Ro + 2f 
8 Ro+1 



(14) 



This relation, which is the central result of this Letter, is 
illustrated in Figure 1. Let us start at the LLL, i.e. at 



Rb 



-1, Rolll = Ro, 



opt V 



-1) 



-0.828. With increas- 



ing Rb, RoQpt(Rb) also increases and reaches the Kep- 
lerian value Ro = -3/4 at Rb = -25/32 = -0.78125. 
At Rb = —1/2 we arrive at solid body rotation, i.e. 
Ro — 0. Interestingly, being connected at Rb = —0.5 
to the branch Ro^p^(Rb) the threshold continues even 
into the positive Ro region corresponding to an outward 
increasing angular frequency. Finally it meets the ULL 
at Ro^pt(— 1) w +4.828 when Rb comes back to —1. 
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1: Dependence of the optimal critical Rossby numbers 
on Rb when Pm = and N — >■ 0. 



Having thus seen that HMRI can easily extend to Ke- 
plerian profiles, we still have to confirm that the shal- 
low B^(t) profiles can indeed be produced by induc- 
tion effects for which some finite value of Rm is still 
necessary. For the sake of illustration, we choose now 
RoKcp — —3/4, and Ha = 30. Figure 2 shows two groups 
of critical curves in the /? — Pm plane. The four curves 
on the right side correspond to SMRI, the curves contin- 
uing into the left part correspond to HMRI. The latter 
ones consist, in general, of two parts, one reaching the in- 
ductionless Pm = area. The connection between them 
typically happens at Rm oc 1. In Figure 3 we show that 
this mechanism is not restricted to n = m = but can 
easily extend to the range of AMRI with higher azimuthal 
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the low Pm-region, as seen in Figs. 2 and 3, may eas- 
ily trigger a quasi-oscillatory behaviour in the non-linear 
regime. Our results encourage experiments on the com- 
bination of MRI and current driven instabilities as they 
are presently planned in the framework of the DRESDYN 
project 01 ■ 
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FIG. 2: SMRI and HMRI for n = 0, Ro = -3/4, Ha = 30. 
Black: Rb = -0.74, Re « 19876, red: Rb = -0.75, Re ^ 
21294, green: Rb = -0.77, Re ki 23935, blue: Rb = -0.91, 
Re 38553, middle of the blue curve Pm = 0.0002, S = 
HaPm^/^ « 0.42, Rm = RePm « 7.71. 
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FIG. 3: Domains of SMRI, HMRI, and AMRI for Ha = 30, 
Ro = -3/4, Re = 4000, Rb = -0.755, and (a) black: n = -1, 
blue: -2, green: -3, red: -4 and (b) red: n = 0, black: 1. 



wavenumbers to, both for small absolute values (Fig. 3a) 
and large absolute values (Fig. 3b) of /3. 

In summary, we have found that the range of appli- 
cability of the inductionless versions of MRI that were 
previously thought to be restricted to Ro < Rolll ~ 
—0.828, can easily extend to Keplerian profiles if only Rm 
is large enough to produce a B^{r) profile that is some- 
what shallower than 1/r. Interestingly, the Ro^pt(Rb) 
curve starting with the ULL further continues to meet 
the ROopt(Rb) branch at the solid body rotation. Since 
this extension of the inductionless forms of MRI circum- 
vents the usual demand S* oc 1, our finding may have 
significant consequences for the working of MRI in the 
colder parts of accretion disks. A detailed investigation 
of the respective roles of 5* and Rm for the onset and the 
saturation mechanism of the instability in different as- 
trophysical problems goes beyond the scope of this letter 
and must be be left for future work. We only note here 
that the sensitive structure of the instability domains in 
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